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(57) Abstract: A system for detecting or determining the nature of a subterranean reservoir. An electromagnetic field is applied 
Q using a dipole antenna transmitter and this is detected using a dipole antenna receiver. The measurements are taken with the antenna 
both in-line and parallel and the difference between the two sets of measurements is exploited. A characteristic difference indicates 
a high resistive layer, which corresponds to a hydrocarbon reservoir. 
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Method and Apparatu s for Determhmipr 
the Nature of Subterranean Reservoirs 

The present invention relates to a method and apparatus for determining the nature of 
5 siibrnarirje and subterranean reservoirs. The invention is particularly suitable for 
determining whether a reservoir, whose approximate geometry and location are 
known, contains hydrocarbons or water, though & can also be applied to detecting 
reservoirs with particular characteristics. 

10 Currently, the most widely used techniques for geological surveying, particularly in 
submarine situations, are seismic methods. These seismic techniques are capable of 
revealing the structure of the subterranean strata with some accuracy. However, 
whereas a sei sm ic survey can reveal the location and shape of a potential reservoir, it 
cannot reveal the nature of the reservoir. 

15 

The solution therefore is to drill a borehole into the reservoir. However, the costs 
involved in drilling ah exploration well tend to be in the region of £25m and since 
the success rate is generally about 1 in 10, this tends to be a very costly exercise. 

20 It is therefore an object of the invention to provide a system for deterrnining, with 
greater certainty, the nature of a subterranean reservoir without the need to sink a 
borehole. 

It has been appreciated by the present applicants that while the seismic properties of 
25 oil-filled strata and water-filled strata do not differ significantly, their 
electromagnetic resistivities (permittivities) do differ. Thus, by using an 
electromagnetic surveying method, these differences can be exploited and the success 
rate in predicting the nature of a reservoir can be increased significantly. This 
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represents potentially an enormous cost saving. 

Consequently, a method and apparatus embodying these principles from the basis of 
the present applicant's co-pending International patent application PCT/GB01/00419. 

5 

This contemplates a method of determining the nature of a subterranean reservoir 
whose approximate geometry and location are known, which comprises: applying a 
time varying electromagnetic field to the strata containing the reservoir, detecting the 
electromagnetic wave field response; seeking in the wave field response, a 
1 0 component representing a refracted wave from hie hydrocarbon layer, and 

determining the content of the reservoir, based on the presence or absence of a wave 
component refracted by die hydrocarbon layer. 

It also contemplates a method for searching for a hydrocarbon containing 
IS subterranean reservoir which comprises: applying a time varying electromagnetic 
field to subterranean strata; detecting the electromagnetic wave field response; 
seeking, in die wave field response, a component representing a refracted wave; and 
. determining the presence and/or nature of any reservoir identified based on the 
presence or absence of a wave component refracted by hydrocarbon layer. 

20 

. It further contemplates an apparatus for determining the nature of a subterranean 
reservoir whose approximate geometry and location are known, or for searching for 
a hydrocarbon containing subterranean reservoir, the apparatus comprising: means 
for applying a time varying electromagnetic field to the strata containing die 
25 reservoir; means for detecting the electromagnetic wave field response; and means 
for seeking, in the wave field response, a component representing a refracted wave, 
thereby enabling the presence and/or nature of a reservoir to be determined. 
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A refracted wave behaves differently, depending on the nature of the stratum in 
which it is propagated. In particular, the propagation losses in hydrocarbon stratum 
are much lower man in a water-bearing stratum while die speed of propagation is 
much higher. Thus, when an oil-bearing reservoir is present, and an EM field is 
5 applied, a strong and rapidly propagated refracted wave can be detected. This may 
therefore indicate the presence of the reservoir or its nature if its presence is already 
known. 

Electromagnetic surveying techniques in themselves are known. However, they are 
10 not widely used in practice. In general, the reservoirs of interest are about 1 km or 
more below the seabed. In order to carry out electromagnetic surveying as a stand 
alone technique in these conditions, with any reasonable degree of resolution, short 
wavelengths are necessary. Unfortunately, such short wavelengths suffer from very 
high attenuation. Long wavelengths do not provide adequate resolution. For these 
IS reasons, seismic techniques are preferred. 

However, while longer wavelengths applied by electromagnetic techniques cannot 
provide sufficient information to provide an accurate indication of the boundaries of 
the various strata, if the geological structure is already known, they can be used to 
20 determine the nature of a particular identified formation, if the possibilities for the 
nature of that formation have significantly differing electromagnetic characteristics. 
The resolution is not particularly important and so longer wavelengths winch do not 
suffer from excessive attenuation can be employed. 

25 The resistivity of seawater is about 0.3 ohm-m and that of the overburden beneath 
the seabed would typically be from 0.3 to 4 ohm-m, for example about 2 ohm-m. 
However, the resistivity of an oil reservoir is likely to be about 20-300 ohm-m. This 
large difference can be exploited using the techniques of the present invention. 
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Typically, the resistivity of a hydrocarbon-bearing formation will be 20 to 300 times 
greater than water-bearing formation. 

Due to the different electromagnetic properties of a gas/oil bearing formation and a 
S water bearing formation, one can expect a reflection and refraction of the transmitted 
field at the boundary of a gas/oil bearing formation. However, the similarity 
between the properties of the overburden and a reservoir containing water means that 
no reflection or refraction is likely to occur. 

10 Thus, an electric dipole transmitter antenna on or close to the sea floor induces 
electromagnetic (EM) fields and currents in the sea water and in the subsurface 
strata. In the sea water, the EM-fields are strongly attenuated due to the high 
conductivity in die saline environment, whereas die subsurface strata with- less 
conductivity potentially can act as a guide for the EM-fields (less attenuation). If the 

IS frequency is low enough (in the order of 1 Hz), die EM-waves are able to penetrate 
deep into the subsurface, and deeply buried geological layers having higher electrical 
resistivity than the overburden (as e.g. a hydrocarbon filled reservoir) will affect the 
EM-waves. Depending on the angle of incidence and state of polarisation, an EM 
wave incident upon a high resistive layer may excite a ducted (guided) wave mode in 

20 the layer. The ducted mode is propagated laterally along the layer and leaks energy 
back to the overburden and receivers positioned on the sea floor. The term 
"refracted" wave in this specification is intended to refer to this wave mode. 

Both theory and laboratory experiments show that the ducted mode is excited only 
25 for an incident wave with transverse magnetic (TM) polarisation (magnetic field 
perpendicular to the plane of incidence) and at angles of incidence close to the 
Brewster angle and die critical angle (die angle of total reflection). For transverse 
electric (TE) polarisation (electric field perpendicular to the plane of incidence) the 
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ducted mode will not be excited. Since the induced current is proportional to die 
electric field, the current will be parallel to the layer interfaces for TE polarisation 
but, for TM polarisation, there is an appreciable current across the layer interfaces. 

5 A horizontal dipole source on the sea floor will generate bom TE and TM waves, but 
by varying die orientation of die receiver antennae, it is possible to vary die 
sensitivity to the two modes of polarisation. It appears that an in-line orientation 
(source and receiver dipoles in-line) is more sensitive to the TM mode of 
polarisation, whereas a parallel orientation (source and receiver dipoles in parallel) is 
10 more sensitive to die TM mode of polarisation. The TM mode is influenced by the 
presence of buried high resistive layers, whereas the TE mode is not By measuring 
with the two antenna configurations and exploiting die difference between the two 
sets of measurements, it is possible to identify deeply buried high resistivity zones, 
i.e. a hydrocarbon reservoir. 

15 

The present invention has arisen from this realisation. 

According to one aspect of the present invention, there is provided, a method of 
determining the nature of a subterranean reservoir which comprises: deploying an 

20 . electric dipole transmitter antenna with its axis generally horizontal; deploying an 
electric dipole receiver antenna in-line with die transmitter; applying an 
electromagnetic (EM) field to die strata containing the reservoir using the 
transmitter; detecting the EM wave field response using the receiver and identifying 
in die response a component representing a refracted wave from die reservoir 

25 according to a first mode; deploying an electric dipole receiver antenna parallel to 
the transmitter; applying an EM field to die strata using the transmitter; detecting die 
EM wave field response using die receiver and identifying in the response a 
component representing a refracted wave from the reservoir according to a second 
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mode; and comparing die first mode retractive wave response with the second mode 
refracted wave response in order to determine the nature of the reservoir. 

According to another aspect of die present invention (here is provided, a method of 
5 searching for a hydrocarbon-containing subterranean reservoir which comprises: 
deploying an electric dipole transmitter antenna with its axis generally horizontal; 
deploying an electric dipole receiver antenna in-line with die transmitter; applying an 
EM field to subterranean strata using the transmitter; detecting the EM wave field 
response using the receiver; seeking in the response a component representing a 
10 refracted wave according to a first mode, caused by a bigh-resistrvity zone; 
deploying an electric dipole receiver antenna parallel to the transmitter; applying an 
EM field to the strata using the transmitter; detecting the EM wave field response 
using the receiver; seeking in die response a component representing a refracted 
wave according to a second mode; and comparing the first mode refractive wave 
15 response with the second mode refractive wave response in order to determine the 
presence and/or nature of any high-resistivity zone. 

_ The first mode may be considered to be a TM mode, and the second mode a TE 
mode. 

20 

Thus, according to the invention, measurements are taken with the transmitter and 
receiver both in-line and parallel, and the two sets of measurements are compared. 
A characteristic difference in values indicates a highly resistive layer located beneath 
highly conductive strata. High resistivity indicates the presence of hydrocarbons and 
25 so the difference in values is a direct hydrocarbon indicator. 



This technique can be used in conjunction with conventional seismic techniques to 
identify hydrocarbon reservoirs. 
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Preferably, the transmitter and/or receiver comprises an array of dipole antennae. 

The technique is applicable in exploring land-based subterranean reservoirs but is 
5 . especially applicable to submarine, in particular sub-sea, subterranean reservoirs. 
Preferably the field is applied using one or more transmitters located on the earth's 
surface, and the detection is carried out by one or more receivers located on the 
earth's surface. In a preferred application, the transmitter's) and/or receivers are 
located on or close to the seabed or the bed of some other area of water. 

10 

In a preferred arrangement, the transmitter and receiver antennae are located on a 
common cable towed behind a vessel. This will result in a fixed offset or a series of 
fixed offsets where several receivers are employed. Preferably, the transmitter 
transmits both modes and may therefore comprise two dipoles arranged mutually at 
IS right angles. Preferably each receiver comprises two dipoles mutually at right 

angles. Preferably one transmitter dipole and one receiver dipole are arranged at 
right angles to the direction of the cable. Alternatively die transmitter and/or 
receivers may escb comprise a single dipole antenna arranged obliquely, eg. at 45° 
to the direction of the cable. With this arrangement the transmitted field is resolved. 

20 

Using mis technique, it is possible to achieve comparable results from the two modes 
as the same signal and offset are used. It will not matter greatly if the transmitter 
drifts in frequency or amplitude. Furthermore, reservoirs can be detected in real 
time. Thus, if the results show a difference in die two modes, this will strongly 
25 indicate the presence of an H/C bearing reservoir and so a more detailed study can 
be made at once. 



Such a system would generally use a single transmission source and several 
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receivers, typically more than ten. The different offsets would be suitable for 
detecting reservoirs at different depths. 

The receivers can be deployed on a single cable or on a series of parallel cables. 
5 There may also be several transmitters. 

In practice, the vessel would normally stop and the cable allowed to sink prior to 
transmission. There would be a transmission at several different frequencies before 
moving to another location. The technique is particularly suitable for edge detection, 
10 and it is a simple matter to select a suitable resolution. However, if the surveying is 
being carried out in an undetermined area, the resistivity in the top layers should be 
mapped, for example with MT methods or by inversion after a reflection study. 

If the offset between the transmitter and receiver is significantly greater than three 
IS times the depth of the reservoir from the seabed (i.e. the thickness of the 

. overburden), it win be appreciated that the attenuation of tite refracted wave will 
often be less than that of direct wave and the reflected wave. The reason for this is 
the fact that the path of the refracted wave will be effectively distance from the 
transmitter down to the reservoir i.e. the thickness of the overburden, phis the offset 
20 along the reservoir, plus the distance from the reservoir up to the receivers i.e. once 
again die thickness of the overburden. 

The polarization of the source transmission will determine how much energy is 
transmitted into the oil-bearing layer in the direction of the receiver. A dipole 
25 antenna is therefore the selected transmitter. In general, it is preferable to adopt a 
dipole with a large effective length. The transmitter dipole may therefore be 100 to 
1000 meters in length and may be towed in two orthogonal directions. The receiver 
dipole optimum length is determined by the thickness of the overburden. 
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The transmitted field may be pulsed, however, a coherent continuous wave with 
stepped frequencies is preferred. It may be transmitted for a significant period of 
time, during which the transmitter should preferably be stationary (although it could 
5 be moving slowly), and the transmission stable. Thus, the field may be transmitted 
for a period of time from 3 seconds to 60 minutes, preferably from 3 to 30 minutes, 
for example about 20 minutes. The receivers may also be arranged to detect a direct 
wave and a wave refracted from the reservoir, and the analysis may include 
extracting phase and amplitude data of the refracted wave from corresponding data 
10 from me direct wave. 

Preferably, the wavelength of the transmission should be in the range 

0.1s £ X <S 5s; 

15 

where X is the wavelength of the transmission through the overburden and s is the 
distance from the seabed to the reservoir. More preferably X is from about 0.5s to 
2s. The transmission frequency may be from 0.01 Hz to 1 kHz, preferably from 1 
to 20 Hz, for example 5 Hz. 

20 

Preferably, the distance between the transmitter and a receiver should be in me range 
0.5X£L^10X; 

25 where X is the wavelength of the transmission through the overburden and L is die 
distance between the transmitter and the first receiver. 

It will be appreciated that the present invention may be used to determine the position, 
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the extent, the nature and the volume of a particular stratum, and may also be used to 
detect changes in these parameters over a period of time. 



The present invention also extends to a method of surveying subterranean measures 
5 which comprises; performing a seismic survey to determine the geological structure of 
a region; and where that survey reveals the presence of a subterranean reservoir, 
subsequently performing a method as described above. 

The invention may be carried into practice in various ways and will now be illustrated 
10 in the following embodiments and reduced scale investigations and simulations. Lathe 
accompanying drawings, 

Figure 1 is a vertical cross-section through a testing tank; 
15 Figure 2 is a plan view of the tank of Figure 1; 

Figure 3 is a plan view of the antennae used in the tank of Figure 1; 
Figure 4 is a side view of the antenna in Figure 3; 

20 

Figures 5 and 6 are respectively a schematic plan view and side view of the testing 
tank set up for measurement; 



Figures 7 is a graph showing calculated and measured values for the transmitted 
25 electric field for a given frequency in the model experiment; 

Figure 8 is a graph showing calculated values tor the electric field in a realistic earth 
model; 



30 



Figure 9 is a schematic side view of a cable layout towed by a vessel; 
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Figure 10 is a plan view corresponding to Figure 9; and 

Figures 1 1 and 12 are views similar to Figure 10 showing two alternative 
5 arrangements. 

The tank 11 shown in Figure 1 and 2 comprises a concrete enclosure 9m long, 6m 
wide and 8m in depth. The tank 1 1 is filled with sea water 12. A diaphragm 13 filled 
with fresh water 14 is located in the tank. The diaphragm 13 is 7.5m long, 4.25m wide 
10 and 0.25m thick and can be located at any desired height in a horizontal orientation 
within the tank 11. 

The conductivity of the sea water 12 was measured to be 5.3 S/m at 14°C and the 
conductivity of the fresh water was measured to be 0.013 S/m. The ratio of the two 
1 5 conductivities is therefore very close to 400. 

The critical frequency f c of a conducting medium, i.e. the frequency at which the 
displacement current is equal to the conduction current, is given by 

20 /c = — — «18 — GHz 

2xSrS» Br 

where Sr is the relative dielectric constant of the medium, and a the conductivity in 
S/m. For water, e, = 80 at the frequencies and temperatures of interest For the two 
25 conductivity values o=* 5.2 S/m and o = 0.013 S/m, ^ = 1.2 GHz and 3 Mhz, 

respectively. Since, in the experiments, the highest frequency is 0.83 MHz, it is a fair 
approximation to neglect the displacement current, even for the fresh water. 



For a non magnetic, conductive medium, the propagation constant v is given by 
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5 

The wavelength X, defined as die distance in which the phase changes In, is given by 
2* /To" 

10 

X m,/in MHz and a in S/m. The skin depth, the distance in which the amplitude 
diminishes by 1/e, is related to the wavelength by 

s= ± = I * 
2x \ 2ffto(T 

15 

For the extremes of the frequency range, for the sea water with o = 5.2 S/m, 

Frequency 30 kHz 830 kHz 

Skin depth 1.27 m 0.24 m 

20 Wavelength 8.01m 1.52 m 

and for the fresh water with a = 0.013 S/m, 

Frequency 30 kHz 830 kHz 

25 Skin depth 25.4 m 4.8 m 

Wavelength 1602 m 30.4 m 
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. Referring now to Figures 3 and 4, two identical electrical dipole antennae, as shown 
were used for the transmitter and receiver. 

Each antenna 15 comprises two square brass plates 16, 15cm square, mounted on an 
5 epoxy substrate 17. Each plate 16 is connected to a co-axial cable 18, which passes 
through an epoxy tube 19 mounted at right angles to the plate 16, to a balun which 
transforms the impedance of the antenna IS from about 2£2 in sea water to about 50& 

The measurement set-up is shown in Figure 5 and 6. An automatic network analyser 
10 (ANA) measures the transmission between the antennae 1 5 as a function of distance 
(offset) and frequency.. The arrangement shown in Figure S shows the antennae IS in 
the parallel orientation. The in-line orientation is achieved by rotating bom antennae 
through 90° in the horizontal plane. 

IS The results of the measurements are shown in Figure 7 together with the 

corresponding theoretical results. The measurements agree well with the theoretical 
results and the figure contains two sets of curves, one with parallel antennae and one 
with the antennae in line. The theoretical results are computed for infmitesirnal dipole 
antennae. The orientations of the antennae and the frequency are shown on the 

20 Figures. 

The parameters of the experiment are scaled relative to possible practical situations. 
To give an idea of orders of magnitude; if the frequency is scaled down by a factor of 
40,000 and the conductivity by.a factor of 10, the dimensions will be scaled up by a. 
25 factor of 632, and the experimental setup would correspond to a low conductivity layer 
of thickness 150 m and conductivity 0.0013 S/m below an overburden of thickness 
300 m and conductivity 0.52 S/m. The corresponding frequency range would be from 
0.75 Hz to 20 Hz, and the length of the antenna nearly 300 m. 

30 The method with the TE and TM-mode have been tested by computer simulations on a 
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simple horizontally layered earth model with electrical parameter values tor typical 
deep water subsurface sediments. The model has an infinite insulating air layer, a 
1150 meter water layer of 031250m, 950 meter overburden of IQm, a 150 meter 
reservoir zone of 50Qm and an infinite underhurden of l£im. Figure 8 illustrates the 
5 amplitude response | E | (electric field) as a function of receiver offset, caused by a 1 
Hz signal. Responses from bom the TM-mode (solid with x's) and TE-mode (dashed 
with +'s) are shown. The amplitudes for the TM-mode are approximately 10 times 
larger at an offset of 5 km. As a reference, the response from a homogenous half- 
space of IQm is shown for both configurations (corresponding to a response from a 
10 water filler reservoir or outside the reservoir area). The TE-mode has the largest 

deviation from its half-space, ie. mis mode is more sensitive to a hydrocarbon layer. 

Figures 9 and 10 show a vessel 31 towing a cable (or streamer) 32 just above the 
seabed 33. The cable 32 carries a transmitter dipole antenna 34 and several receiver 
15 dipoles 35, only four of which are shown. The depth of water might be of the order of 
1 000m, the offset between the transmitter 34 and the nearest receiver 35 might be 
about 2000m and the receivers might be about 100m apart The transmitter 34 is 
controlled from the vessel 3 1 via the cable 32 and the responses detected by the 
receivers 35 are relayed back to the vessel' 3 1 in real time, again via the cable 32. 

20 

Figure 10 shows an arrangement in which the vessel 31 tows three cables 41, 42, 43, 
each carrying a series of receivers 45, 46, 47. The spacing of the three cables 41, 42, 
43 is achieved by means of a spar 44. 

25 In the arrangement shown in Figure 1 1, the transmitter 48 is in the form of two dipole 
antennae, one parallel to the cable 42 and one at right angles. 

The arrangement shown in Figure 12 is similar to Figure 1 1, but in this case, the 
transmitter 5 1 is a single dipole antenna arranged at 45° to the cable 42. 

30 
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Claims 

1. A method of determining the nature of a subterranean reservoir which 
comprises: deploying an electric dipole transmitter antenna with its axis generally 

5 horizontal; deploying an electric dipole receiver antenna in-line with the transmitter; 
applying an electromagnetic (EM) field to the strata containing the reservoir using 
the transmitter; detecting the EM wave field response using the receiver and 
identifying in the response a component representing a refracted wave from the 
reservoir according to a first mode; deploying an electric dipole receiver antenna 

10 parallel to the transmitter; applying an EM field to the strata using the transmitter; 
detecting the EM wave field response using the receiver and identifying in the 
response a component representing a refracted wave from the reservoir according to 
a second mode; and comparing the first mode refractive wave response with the 
second mode refracted wave response in order to determine the nature of the 

IS reservoir. 

2. A method of searching for a hydrocarborHX>nraining subterranean reservoir 
which comprises: deploying an electric dipole transmitter antenna with its axis 
generally horizontal; deploying an electric dipole receiver antenna in-line with the 

20 transmitter; applying an EM field to subterranean strata using the transmitter; 
detecting the EM wave field response using the receiver, seeking in the response a 
component representing a refracted wave according to a first mode, caused by a 
high-resistivity zone; deploying an electric dipole receiver antenna parallel to the 
transmitter; applying an EM field to the strata using the transmitter; detecting the 

25 EM wave field response using the receiver; seeking in the response a component 
representing a. refracted wave according to a second mode; and comparing the first 
mode refractive wave response with the second mode refractive wave response in 
order to determine the presence and/or nature of any high-resistivity zone. 
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3. A method as claimed in Claim 1 or Claim 2, characterised in that the first 
mode is a TM mode of polarisation and/or the second mode is a TE mode of 
polarisation. 

5 

4. A method as claimed in any preceding Claim, characterised in that the 
transmitter and/or receiver comprises an array of dipole antennae. 

5. A method as claimed in any preceding Claim, characterised in that the 
10 transmitter and/or receiver is located on or close to the seabed or the bed of some 

other area of water. 

6. A method as claimed in any preceding Claim, characterised in mat the 
transmitter and receivers are located on a common cable arranged to be towed 

IS behind a vessel. 



7. A method as claimed in any preceding Claim, characterised in that die 
transmitter comprises two dipole antennae arranged mutually at right angles. 

20 8. A method as claimed in any preceding Claim, characterised in that each 
receiver comprises two dipole antennae arranged mutually at right angles. 

9. A method as claimed in Claim 6 characterised in that the transmitter and/or 
receiver each comprise a single dipole antenna arranged obliquely to the direction of 

25 the cable. 

10. A method as claimed in any preceding Claim, characterised in that the 
frequency of the EM field is continuously varied over die transmission period. 
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11. A method as claimed in any preceding Claim, characterised in mat the field is 
transmitted for a period of time for 3 seconds to 60 minutes. 

5 12. A method as claimed in Claim 11, characterised in that the transmission time 
is from 3 to 30 minutes. 

13. A method as claimed in airy preceding Claim, characterised in that the 
wavelength of the transmission is given by the formula 

0.1s < X <10s; 

wherein X is the wavelength of the transmission through the overburden and s is the 
distance from the seabed to die reservoir. 

14. A method as claimed in any preceding Claim, characterised in that distance 
between die transmitter and a receiver is given by the formula 

0.5<X<10s; 

20 

where X is the wavelength of the transmission through the overburden and 1 is die 
distance between the transmitter and the receiver. 

15. A method as claimed in any of Claims 10 to 14, characterised in that die 
25 transmission frequency is from 0.01 Hz to 1 kHz. 

16. A method as claimed in Claim IS, characterised in that the transmission 
frequency is from 1 to 20 Hz. 



10 



15 
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17. A method as claimed in any preceding Claim, characterised in mat it includes 
suppressing the direct wave and/or any other known wave contribution that may 
disturb the measurements, thereby reducing the required dynamic range of the 
receiver and increasing the resolution of the refracted wave. 

18. A method of surveying subterranean measures which comprises: performing 
a seismic survey to determine the geological structure of a region and where that 
survey reveals the presence of a subterranean reservoir, subsequently performing a 
method as claimed in any preceding Claim. 
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